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Abstract 

We study the coupled rotation-vibration levels of a hydrogen molecule in a confining potential 
with cylindrical symmetry. We include the coupling between rotations and translations and show 
how this interaction is essential to obtain the correct degeneracies of the energy level scheme. We 
applied our formalism to study the dynamics of H2 molecules inside a "smooth" carbon nanotube 
as a function of tube radius. The results are obtained both by numerical solution of the (2 J + 1)- 
component radial Schrodinger equation and by developing an effective Hamiltonian to describe 
the splitting of a manifold of states of fixed angular momentum J and number of phonons, N. 
For nanotube radius smaller than ~ 3.5 A, the confining potential has a parabolic shape and the 
results can be understood in terms of a simple toy model. For larger radius, the potential has the 
"Mexican hat" shape and therefore the H2 molecule is off-centered, yielding radial and tangential 
translational dynamics in addition to rotational dynamics of H2 molecule which we also describe 
by a simple model. Finally, we make several predictions for the the neutron scattering observation 
of various transitions between these levels. 

PACS numbers: 78.70.Nx,34.50.Ez,82.80.Gk,71.20.Tx,36.20.Ng,63.22.+m 



I. INTRODUCTION 



The study of quantum dynamics of hydrogen molecules in confined geometries has recentl 
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developed into an active field both experimentally and theoretically 

Ifl; lllll^l due to potential use as catalysts, molecular sieves, and storage media. In the case 
of fullerenes and nanotubes, such trapping may yield new exotic quantum systems due to 
zero and one dimensionality of the absorption sites, respectively. Thus, understanding the 
structural and dynamical aspects of trapping in confining geometries is of both fundamental 
and practical importance. 

The theory of molecular rotation in solids has a long history dating back to the early work 



of Pauling 



1 31 ] . Devonshire [1J], and Cundy|l5j]. They introduced the concept of the crystal 



field potential, V(Q), where Q specifies the orientation of the molecule, to solve for the 
energy levels of the hindered rigid- rotor. This traditional approach assumes that the center 
of mass (CM) of the trapped molecules are fixed and therefore does not take into account 

n n 

the rotation- vibration (RV) coupling. However recent studies p, have indicated that 
vibrational levels of H 2 trapped in the octahedral sites of C 60 , for example, are significantly 
perturbed by RV coupling and in a previous paper [5] (I) we have shown that this coupling 
has to be included in a symmetry analysis of the energy level degeneracies. Interestingly, 
to date there is a little done to treat CM dynamics and RV coupling. Most of the studies 
are based on the approximation where an effective orientational crystal field potential is 
obtained after the potential is averaged over the zero-point translational motions of the H2 
molecule jll. ll^ . 

Recently (in I) we have presented a detailed theory of coupled RV dynamics of H 2 molecule 
trapped in a zero dimensional cavity with various symmetries. Here we present a similar 
study to analyze the combined rotational and translation states of hydrogen molecules con- 
fined in a one dimensional potential. This problem is closely related to the experimental 
situation where hydrogen molecules are absorbed into carbon nanotube ropes. Figure 1 
shows schematically various types of absorption sites for H2 molecule. Several neutron and 
Raman scattering experiments have been carried out to characterize the binding ene rgie s 
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and rotational barriers for H2 at these sites with conflicting results^, 0, U, 
One of the motivations of the present work is to provide a detailed description of the RV 
dynamics of H2 molecules at these different absorption sites and discuss the consequences 
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for inelastic neutron scattering experiments. In the present paper, we focus our attention 
on the general formalism and discuss only the case where a single H2 is confined inside a 
single nanotube. Extension of this work to the interstitial and external sites and to cases 
where H 2 molecules interact with one another will be presented elsewhere. 

Briefly this paper is organized as follows. In the next section, we discuss the potential 
model for hydrogen and nanotube interactions and validate several approximations, such as 
assuming a smooth tube, used in our formalism. In Sec. Ill we present our formalism to 
treat the coupled rotational and translational motion of H 2 molecule confined in a smooth 
nanotube. We show that the problem can be mapped into a (2 J + l)-component radial 
Schrodinger equation which can be solved numerically. In Sec. IV we discuss the dynamics 
of a hydrogen molecule when the confining potential has parabolic shape (which occurs for 
a small-radius nanotube). We interpret the exact numerical results in terms of a simple 
analytical toy model. In Sec. V we discuss the case where the confining potential has a 
Mexican-hat shape. For this case (which occurs for large radius nanotubes such as (10,10)) 
the equilibrium position of the CM of the H2 molecule is off-center and it performs radial 
and tangential translational oscillation in combination with its rotational dynamics. In 
this section, we also present several perturbation results which help to interpret the exact 
numerical results. In Sec. VI we discuss the experimental observation of various transitions 
via inelastic neutron scattering measurements. Our conclusions are summarized in Sec. VII. 

II. POTENTIAL MODEL 

We model the intermolecular potential for H2 trapped in a carbon nanotube as a sum of 
atom-atom potentials 



The dependence of the potential on the position (r) and orientation (Q) of H2 molecule is 
through the inter-atomic distances r^. All the results reported in this paper are obtained 
from the same WS77 potential, 1 6] -A/r 6 + B exp(-CV), (where A = 5.94 eVA 6 , B = 678.2 



the WS77 potential gave the best fit to the energy spectrum of H 2 in solid C 50 . 

For simplicity we will restrict this formulation to the idealized case when the hydrogen 




i,H j,C 



(1) 




Compared to other commonly used potentials, 
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molecule is confined by a so-called "smooth" nanotube. By this we mean that the poten- 
tial produced by the nanotube has cylindrical symmetry and is invariant with respect to 
translations along its axis of symmetry. In some of our numerical work we will study "real" 
nanotubes which do not posses the high symmetry of "smooth" nanotubes. 

It is instructive to look at various potentials for an orientationally averaged hydrogen 
molecule (i.e. para hydrogen with J = 0) when H 2 is inside and outside a single nanotube. 
Figure 2 indicates two different types of confining potential depending on the nanotube 
radius. Figure 2a shows that for small nanotubes such as (9,0), the potential minimum 
occurs at the center of the tube and therefore the potential has a parabolic shape. However 
for larger nanotubes such as (10,10) nanotube, the minimum is off-centered and therefore 
the potential has a Mexican-hat shape. Because of this off-centering the dynamics of the 
H 2 molecule is a quite interesting and rich one as we discuss in detail below. The right 
panels in Fig. 2 shows the potential when the H 2 molecule is outside the nanotubes. The 
outside binding energy does not depend on the tube radius strongly and is about 30 meV. 
The horizontal lines indicates the radial phonon energy levels, indicating that at least a few 
bound states can occur even for a hydrogen molecule outside a single nanotube. The solid 
and dashed lines in Fig. 2 shows the results with and without the smooth tube approxima- 
tion, respectively. Since these two curves are very close to one another, the smooth tube 
approximation will not cause significant error in our theory. 

In order to develop some intuition about the orientational potential for a hydrogen 
molecule in a nanotube, in Fig. 3a we show the radial potential for three different orienta- 
tions of H 2 molecule inside a (10,10) nanotube. We note that the radius of the (10,10) tube 
is large enough that the parallel to the axis (p) and tangential (t) orientations (as depicted 
in the inset to Fig. 3a) give almost the same energy. However, the radial (r) orientation of a 
H 2 molecule has a minimum energy which is about 8 meV higher in energy than that of the 
other two orientations. We also point out that the position of the CM of the H 2 molecule 
for the minimum potential energy changes about 0.2 A depending on the orientation of 
the H 2 molecule. This is a clear indication that the orientational and vibrational motion of 
hydrogen molecule are significantly coupled. 

Figure 3b shows the minimum potential energies, E r , E p , and E t for respective radial, 
parallel, and tangential orientations of an H 2 molecule inside various nanotubes. It is clear 
that for nanotube radius around 3 A, the orientational dependence of the potential is of the 
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order of 30 meV and therefore is comparable to the energy separation of 80 meV or more 
between energy levels corresponding to different rotational quantum number J's. Accord- 
ingly, for tubes whose radius is less than that of a (9,0) tube, J may not be considered to be 
a good quantum number. For large nanotubes, such as (10,10), the orientational potential 
is of the order of 8 meV and does not change much with larger radius tubes (i.e. close to 
the graphite limit). In the present paper, we present our formalism using (9,0) and (10,10) 
nanotubes which represent the two potential regimes; namely the parabolic and Mexican-hat 
potentials, but for both tubes J is considered to be a good quantum number. 



III. FORMULATION 



The hydrogen molecule is unique in that its moment of inertia is small enough that the 
rotational kinetic energy often dominates the orientational potential in which the molecule 
is placed. Under these circumstances the rotational quantum number J is nearly a good 
quantum number and the effect of the orientational potential is to reduce the degeneracy 
of the 2 J + 1 sub-states of a given J. (The generalization of the formulation we present 



below to the case when J is not a good quantum number will be presented elsewhere [171]). 
In the present case any eigenfunction describing the orientational and translation state of 
the molecule can be written in the form 

*(r,0) = ^(r,0 r ;fi)e^= £ $^ fc (r, <f> r )Yf (Q)e ikz , (2) 

M=-J 

where r, z, and <p r are the cylindrical coordinates of the center-of-mass of the hydrogen 
molecule, Q denotes its molecular orientation specified by angles 6 and <f), and Y"j w (f^) is 
a spherical harmonic. We will refer to ip(r,(f) r ;Q) as the cylindrical RV wavefunction. For 
economy of notation we henceforth omit the superscripts J and k. Because $m is allowed 
to depend arbitrarily on r and (f) r , this wave function takes into account the most general 
interaction between rotations and translations subject to the constraint that J is a good 
quantum number. In this notation, the Hamiltonian for a single hydrogen molecule (of mass 
m) with J and k fixed is written as 

h 2 

H = -—V 2 + BJ(J+l) + V(r,n), (3) 

where m is the mass of an H2 molecule and V(r, Q) is the orientational potential which, as 
indicated, also depends on position. For a smooth nanotube we may write the orientational 



potential energy as 



V(v, fl) = Vo{r) + £ ^V, <f> r )Y™(e, </>) , (4) 



where the sum is over L > 0. Because the hydrogen molecule is centrosymmetric, only 
terms with L even appear in the potential. Also, because a smooth nanotube has a mirror 
plane perpendicular to the axis of symmetry, only terms with even M appear. Furthermore, 
because a global rotation of the molecule (i. e. incrementing both <fi and r by the same 
amount) is a symmetry of the system, we may write 

V{t, n) = Vo{r) + £ v™(r)YM(9, 0)e^-« , (5) 

L,M 

where vff is a function only of r and % M = f^*. In addition, vff(r = 0) vanishes for 
M 0. There is also a mirror plane containing the long axis of the tube which implies 
that the potential should be an even function of (0 — (fi r ). This implies that vff(r) is a 
real- valued function. This function may be evaluated by integrating the potential at a fixed 
center-of-mass position over all orientations: 

v ^(r) = e iM ^ J dQY^(9, <f>)*V(r, SI) . (6) 

Contrary to appearance, v^f(r) does not depend on <p r because V(r, Q) is a function of 

- <Pr)- 

The Schrodinger equation for ip(r, <ft r ; Q) is 

= # a ty (a) (r,^;n) , (7) 

where = £ (a) - 5J(J + 1) - h 2 k 2 /(2m). For given values of J and fc, this equation 
generates a spectrum of eigenvectors ip^ a \r, <p r \ VL) with associated eigenvalues E^ a \ for a = 
0, 1, . . .. Substituting Eq. (J2J) into the Schrodinger equation we rewrite it in the form 



h 2 d 2 Id 1 d 2 \ .[r a) , 

"2^ [d^ + rd-r + ^M +Vo(r) r ,,(r,W 



L,M' V / 

= ^ a )$£ ) (r,0 r ) , M = — J, — J + 1 ... J — 1, J • (8) 
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We see that we may write a solution to this set of equations in the form 



(9) 



where P is a quantum number whose significance we will discuss shortly. Thus we have 



For each P index we have a Schrodinger equation for the (2J+ l)-component wave function 
which is of the form [/_jp(r), /ij+i p( r ) ; • • • fjp( r )]- Then the cylindrical RV wavefunction 
is 



= e-^Y,fMAr)Yrme~ m ^ = e iP +* £ f$ P (r)Yj M (9)e iM ^ (11) 



It is important to keep in mind that vff(r) vanishes for odd M. As a consequence, in the 
(2 J + l)-component wavefunction there is no mixing between even and odd values of M. 
For J = 1 wavefunctions one will have "even" wavefunctions in which the sum over M in 
Eq. (fTT|) reduces to the single term for M = and "odd" wavefunctions in which the sum 
over M in Eq. (|llj) includes only M — ±1. 

The quantum number P indicates that this wavefunction transforms like e~ tP< ^ when 
the position and orientation of the molecule are simultaneously rotated about the axis of 
symmetry. Under this global rotation the quantity 



is invariant, so the total wavefunction transforms as stated. When vff(r) is independent 
of r, then, since vff(r = 0) must vanish, we have that vff(r) = for M / and, in 
Eq (fTUj) there is no coupling between /j^'s for different values of M. Thus, in this case 
there is no dynamical interaction between the orientational coordinate and the center-of- 
mass coordinate and the wave function can be chosen so that Jm,p is only nonzero for a 
single value of M. Then one has the usual separation of variables so that the orientational 
wave function is proportional to e tM ^ and the translation wave function is proportional to 
e -»(M +P)<t> r _ Here we have accomplished a similar separation of coordinates when vff(r) is 
allowed to depend on r. Now the result is not a scalar radial equation, but rather a radial 




M 



M 



(12) 
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equation for a (2J + l)-component wave function. That is the result embodied in Eq. (fTUj). 
where we have one such (2 J + l)-component radial Schrodinger equation for each value of 
P. 



IV. QUASI-HARMONIC POTENTIAL 

In this section we discuss the case exemplified by a H 2 molecule inside a (9, 0) nanotube, 
for which the minimum of the potential Vo(r) occurs for r = 0, in which case we will introduce 



A. No Interactions between Rotations and Translations 

1. Harmonic Potential 

Here we discuss the eigenvalues and eigenfunctions of the two-dimensional isotropic har- 
monic oscillator, to emphasize the relation between the above formulation in terms of cylin- 
drical coordinates and that in terms of Cartesian coordinates. 

For an isotropic and harmonic potential we expect the eigenvalues to be 



Note that the nth level (with energy nhiS) is ra-fold degenerate, because in Cartesian no- 
tation, if, say n = 4, we have wave functions (3,0), (2,1), (1,2), and (0,3), where (n, m) 
is a wave function with n excitations in the x coordinate and m excitations in the y coor- 
dinate. This degeneracy reflects the U2 symmetry corresponding to the invariance of the 
Hamiltonian with respect to a transformation of the form 



where a' x and a' y create phonons in the x and y-coordinates, respectively and U is a two- 
dimensional unitary matrix. This transformation is essentially the same as a four dimen- 
sional rotational symmetry in the space of the momenta p x -, Py, and coordinates x and y. 
Since the kinetic energy is quadratic in the momenta, spherical symmetry in this space only 
holds if the potential is harmonic. 



a toy model with the isotropic harmonic potential Vo(r) 





(13) 




(14) 
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In cylindrical coordinates the eigenfunctions can be written as 
where i/)%(r) satisfies the radial equation, 



2m 



d 2 r»(r) l#«(r) /x 2 



dr 2 



+ - 

r 



dr 



(15) 

(16) 
in order of 



Here the family of solutions for a given value of /x are labeled a = 0, 1,2, 
increasing energy and for the isotropic and harmonic potential we have 

El = (/x + 1 + 2a)Uuj . (17) 

So from the radial equation for /x — we have eigenvalues Hcu, Shoo, 5hu, etc. The fact that 
we have a seemingly accidental degeneracy between different representations (i. e. between 
different values of /x) is the result of the U2 symmetry of the Hamiltonian mentioned above. 
A consequence of this symmetry is that for a harmonic potential the total energy depends 
only on the total number of phonon excitations. This symmetry is distinct from the circular 
symmetry in x-y space. 



2. Anharmonic Potential 



The U2 symmetry is broken by anharmonic terms which then take us into the generic case 
of a particle in a circularly symmetric potential which is not harmonic. Accordingly, we now 
consider the effect of adding an anharmonic perturbation of the form 7r 4 to the harmonic 
potential. For illustrative purposes, we treat this anharmonic perturbation within first-order 
perturbation theory. Our results are characteristic of the generic case, for which different 
values of m give rise to distinct eigenvalues. In this case, the n-fold degenerate manifold 
which has energy nhu for the harmonic potential is split into doublets (corresponding to 
the degeneracy between +m and — m) and, if n is odd, a singlet from m — 0. Our explicit 
results are given in Table I. These results are generic in the sense that addition of further 
anharmonic terms will not further change the degeneracies. 



B. Toy Model of Translation- Rotation Coupling 



In this section we explore the consequences of allowing coupling between rotations and 
translations. Since we now restrict attention to the manifold of (J — 1), we need keep only 
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terms with L = 2 in Eq. ((SJ). Thus, as a toy model, we set 

V(r, O) = l -kr 2 - ^-a{3 cos 2 9-1)- ^(3r 2 sin 2 9 cos(20 - 20 r ) , (18) 

where a, (3, and k are constants and the factor — ~ is included so that the matrix elements are 
numerically simple. Also we take the dependence on r to be quadratic to facilitate calculation 



of the matrix elements. In the language of Eq. © this model has ^(r) = — q;v207T and 
v f 2 ( r ) = — (3r 2 J lOn/ 3. We are going to consider the effect of this Hamiltonian within 



the manifold of (J = 1) states. Using this toy model we can illustrate how the rotation- 
translation affects the symmetry of the energy levels. Within the (J = 1) manifold we may 
use operator equivalents to write 

V(r,0) = hr 2 + a(3J 2 -2)+(3[(J 2 -J 2 )(x 2 -y 2 ) + 2(J x J y + J y J x )xy] 

= \kr 2 + a(3J 2 - 2) + 1 -(3{J 2 + + J 2 )(x 2 - y 2 ) - i(3{J 2 + - J 2 _)xy . (19) 

For illustrative purposes we will assume that a and j3a 2 are small compared to the phonon 
energy huo. Here a 2 = (x 2 ) = (y 2 ), where the averages are taken in the ground state. In 
that case, in addition to the quantum number P, the total number of phonons, N, is a good 
quantum number. However, we emphasize that in our numerical work 3] , we do not make 
this approximation. Equation (fTU|) assumes that J is a good quantum number but mixes 
states with different numbers of phonon excitations. We should mention that the toy model 
assumes that the molecule has minimal potential energy when it on the axis of the tube. For 
small [e. g. (9,0) tubes] this assumption is justified. For larger tubes, the minimal potential 
energy occurs for a nonzero value of r and the molecule is dominantly off center. We will 
later treat that case using a different model. 

C. Results of the Toy Model 

We now discuss the results of the toy model assuming that the number of phonons is a 
good quantum number. We note that all the energy expressions given below are with respect 
to BJ(J+ 1) with J = 1. 
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1. Zero Phonon Manifold 



We first consider the manifold of the states having J = 1 with zero phonons. One finds 
that the energy is diagonal in J z with 



so that (if a is positive) one has the singlet J z = state lower than the doublet J z — ±1 
states by an energy separation of 3a. One may visualize this as the energy difference between 
a state for which the molecule is in the phonon ground state and is oriented parallel to the 
axis and the two states when the molecule is in the phonon ground state and is oriented 
transversely to the axis. For later use we tabulate these wavefunctions in Table |TT| 

2. One-Phonon Manifold without Rotation- Translation Coupling 

If we set P = in the toy model of Eq. (fTUj) . then essentially we have independent 
oscillation of molecules which have fixed orientation. Then if n x and n y are the vibrational 
quantum numbers, we see that in the one-phonon manifold (n x + n y = 1) we have 



so that the lowest energy state (if a > 0) is doubly degenerate and the excited state is 
four-fold degenerate, as is shown in Fig. 0] 

3. One-Phonon Manifold with Rotation- Translation Coupling 

The unphysical aspect of the energy level scheme we just found for the one-phonon 
manifold is that it does not take into account that the molecular orientation ought to be 
correlated with the translational motion. If the molecule translates near the wall, then the 
molecule should preferentially be parallel to the wall. This means that the orientation of 
the molecule has to be correlated with the translational motion. This effect will be greater 
the more strongly the wall potential affects the motion of the molecule. 

In terms of number operators n x and n y which are the number of phonon excitations in 
the x-direction and ^-direction, respectively and a x and at which are creation operators for 



E(J Z ) = a(3J 2 z - 2) , 



(20) 



E(n x , n y , J z ) = 2hu + a(3J| — 2) , 



(21) 
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these excitations, we may write the Hamiltonian for the one-phonon ( J = 1) manifold as 



H = huo(n x + n y + 1) + ct(3J 2 2 - 2) 



+ o 2 f3 



(J+ + /^(n-c - - i(J+ - J^)(aJ,a 3/ + a^a x ) 



(22) 



This gives the energy level scheme shown in the rightmost panel of Fig. 0] The wavefunc- 
tions are given in Table 11111 and we discuss them now. First of all, in a classical picture, we 
would argue that the molecule can oscillate equivalently in each of the two coordinate direc- 
tions transverse to the cylinder. In each of these two cases the molecule can assume three 
inequivalent orientations because the directions a) along the axis of the tube, b) parallel to 
the directions of spatial oscillation, and c) transverse to the direction of spatial oscillation 
are all inequivalent to one another. This argument predicts that the six states form three 
doubly degenerate energy levels. Quantum mechanically the situation is different. In Fig. 0] 
we show the energy levels when no dynamical mixing between rotations and translations is 
allowed, i. e. for (3 = 0. In this limit the six states form a degenerate doublet and a degen- 
erate quartet. When translation-rotation mixing is allowed, i. e. for (3 ^ 0, we now have 
the generic case of two doublets and two singlets, as shown in Fig. |3J The wavefunctions are 
shown in Fig. (It is interesting to note that it is not obvious that the wave functions for 
P = +2 and P = —2 are related by symmetry.) 



4- Two-Phonon Manifold with Rotation- Translation Coupling 

Actually, because the dependence on r of the matrix elements in Eq. (fTT?|) was taken to 
be either constant or proportional to r 2 , the representation of Eq. (|2*2*|) is valid within any 
manifold of fixed total number of phonons and (J = 1). The removal of the degeneracy in 
the energy level scheme of the two-phonon (i.e. n x + n y = 2) and (J = 1) manifold according 
to the Hamiltonian of Eq. (|22|) is shown in Fig. |HJ The eigenf unctions and eigenvalues for 
this manifold including anharmonicity are listed in Table HV1 

The results obtained by numerically solving the eigenvalue problem of Eq. (fTUJ) using 
WS77 potential are given in the last column of Table IV[l8]. To understand the meaning 
of this spectrum, we relate these results to those of the toy model when the parameters 
of the toy model are suitably chosen. For a good fit we allow the constants a and (3 to 
depend on the total number of phonons, N. (This dependence reflects the fact that the 
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dependence of the parameters of the toy model on r is arbitrary and unrealistic.) In this 
simple model we also include the anharmonic term 7?" 4 which we treat within first order 
perturbation theory. We determine the best parameters for the toy model by making a least 
squares fit of the numerically determined energy levels to those of the toy model and these 
parameters as well as the results of this fit are given in Table. The fact that depends 
on N indicates that we should probably replace a by or 2 . Also the fact that the splitting 
of the two phonon manifold is not perfectly reproduced by the toy model indicates that the 
anharmonicity energy is not simply proportional to r 4 . Nevertheless the close agreement 
between our numerical results and those of the toy model indicates that this model provides 
a useful simple picture of translation-rotation coupling. 

5. Summary 

We can summarize the systematics of the rotation-translation spectrum of the toy model 
we have introduced. We first consider the harmonic 7 = case and then discuss the effect 
of introducing anharmonicity. In the iV-phonon sector the harmonic phonon wavefunctions 
give rise to states proportional to (x + iy) N . Combining these with a J z = 1 state gives a 
unique P = N + 1 state. This state will be degenerate with the similar P = —N — 1 state. 
In the toy model these states have energy (N + l)tkv + a. Adding anharmonicity shifts the 
energy of these two states, but their degeneracy is generic. 

Harmonic phonon states which transform like (x + iy) N ~ 2k will uniquely combine with 
J z = states to form states for which P = (N — 2k) and which have energy (N + l)tkv — 2a. 
In analogy with Fig. 4, anharmonicity splits these states into doublets of +P and — P and, 
if P is even, a singlet for P = 0. 

The rotation-translation coupling (proportional to f3) influences the states with P = 
N—l, P = N — 3, etc. For positive P one has two eigenstates made from linear combinations 
of states of the form <pi = (x + iy) p+1 \J z = —1) and </> 2 = (2 + iy) p ~ l \ J z = +1). Since the 
rotation translation coupling interaction proportional to xy(J+ — J 2 ) has matrix elements 
between these two states, the eigenstates <pi ± <p2 will be split by an amount proportional to 
(3a 2 and this splitting will be modified by anharmonicity. Obviously, this scenario indicates 
that one can not understand the degeneracies of the states of a hydrogen molecule in confined 
geometry without considering the effect of rotation-translation coupling. 
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V. MEXICAN HAT POTENTIAL 



Here we discuss the case when the minimum of the potential Vo(r) occurs for non-zero r 
as happens for H 2 molecules inside 10 x 10 tubes or for H 2 molecule in a bound state outside 
any tube. We start from Eq. (|10|). To see what this equation yields, we first consider 
its solutions for a (J = 0) molecule. We have solved the eigenvalue problem of Eq. ()10|) 
numerically on a mesh of points for a 10 x 10 tube[18]. The results shown in Fig.[7^ indicate 
two different regimes for the dependence of the energy levels on the quantum number P. For 
the low-lying energy states it is quadratic and then gradually becomes linear as the energy 
of the states increase. 

It is possible to understand the quadratic behavior of energy levels versus the quantum 
number P based on a simple idealized model. Assuming Vo(r) can be replaced by a harmonic 
oscillator potential, Eq. (fTUj) becomes essentially 

h 2 d 2 h 2 p 2 



+ n 2 + + ^(r - r 



,2 



f { p\r) = EP ff\r) . (23) 



2m dr 2 2mr 2 2 

In writing this result we dropped the term linear in the derivative. This term does not 
contribute to the energy in first-order perturbation theory. When we treat the term in P 2 
perturbatively, this equation leads to a harmonic oscillator spectrum with 

<» = B „ + ( J V + i)^ + ^/-L\, (24) 



where u = ^k/m and (X) here indicates an average of X over the radial wavefunction. 

The gray solid lines in Fig. shows the results based on this model and the points are 
from the numerical exact results, indicating that our idealized model successfully describes 
the low-lying energy spectrum. Here hu is of order 14 meV and and the quantum of tangen- 
tial kinetic energy (h 2 /(2mr 2 )) is about 0.1 meV. Curiously, this spectrum is reminiscent of 
the vibration- rotation spectrum of a diatomic molecule . Finally as we go away from the 
ground state, the simple model is not enough to explain the observed behavior. We note 
that the spacing between the energy levels is not constant (probably due to an anharmonic 
contribution to the potential) and the dependence on the P becomes almost linear. 

We next discuss the solution of Eq. (fTUJ) for a ( J = 1) H 2 molecule. Fig. 0d shows the 
results obtained numerically, ^j] . We will interprate numerical results using a simple model 
which includes the translation-rotation coupling (as embodied by the vff's). We expect the 
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radial wavefunctions to be Gaussians centered about r = tq. Indeed in the terms containing 
vff(r), we will make the replacement 



<(r)-«(r)> . 



(25) 



(To a good approximation these averages can be calculated from the ( J = 0) wavefunction.) 
In what follows we set 



(v%\r)) = w M 

for M = and M = 2. Then, if we set r = tq + x, Eq. (jlOj) may be written as 



(26) 



h 2 d 2 



2m dr 2 
1 



h 2 P 2 „ 1,, v 2 
+ W +Eo+ 2 fc(r - ro) 



^7T 



^C(121;M-M',M')^/, 



(«) 

M-M'P 



(a) 
PJM,P 



E a P f, 



(27) 



where the Clebsch-Gordan coefficients assume the values C(121; L, 0) = (3L 2 — 2)/vT0 and 
(7(121; 



-1,2) = y3/5. For each value of (the number of radial phonons) and P (the 
number of tangential excitations) the Hamiltonian is the following three dimensional matrix 
(where the rows correspond to J z — — 1, J z — 0, and J z — +1, in that order): 



W(JV,P) = [E + iV^]X + 



15 + A(P 



5 








P 




(28) 



where X is the unit matrix, A = (h 2 /(2mr 2 )) , 5 = — 3u>o/\/207r, and £> = — W2y3/107r. For 
fixed values of N and P we have the three energy eigenvalues 



E (o) 

£(±) 



P + Nhuo --5 + AP 2 
3 



P + Nhu + -5 + A(P 2 + 1) ± V4A 2 P 2 + P 2 
3 



(29) 



In Fig. [TJd we show the spectrum of a (J = 1) molecule obtained numerically from 
Eq. (|Tnj) as a function of P. Our numerical results indicate the phonon number N is a good 
quantum number and can so be identified only for N < 3. Accordingly, we limit our detailed 
interpretation in terms of the model of Eq. (f2*Kj) to = 0. For each value of P there are 
three energy eigenvalues, two of which are close in energy. These corresponds to the case 
where the H 2 molecule is oriented parallel to the tube surface (i.e. t and p orientations in 
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Fig. 3). The third energy corresponds to the orientation perpendicular to the tube surface 

(i.e. r radial orientation shown in Fig. 3). This orientation has an energy about 3 meV near 

than that of the other two orientations, and is comparable to 2.6 meV observed for H 2 on 
I I 

graphite [12| . Fig. [7jp shows also the result of the simple model of Eq. (|28|) . The values of A, 
B, and S used to get a good fit are given in the figure. The value of A (0.092 meV) is not very 
different from the value h 2 / (2mrp) = 0.087 one gets from the value of r = tq = 3.46 A at 
the minimum of the potential. Thus the numerical results are easily understood in terms of 
our simple toy model. 



VI. EXPERIMENTAL OBSERVATION OF THE ENERGY SPECTRUM 



Here we make some remarks concerning the observation of these modes via inelastic 
neutron scattering. Specifically we consider the energy loss spectrum in the neutron time- 
of-flight spectrum. (This technique has been used to probe local excitation of H 2 molecules 
in the octahedral sites of C6o-^|) We start be recalling the results for the cross section for 
inelastic neutron scattering of H 2 molecules. When the very small coherent (i. e. nuclear 
spin independent) scattering is neglected, the result is 

d 2 a k' 



dVtdE k 



NxS^ + NxS^o + N(l - x)S ^ 1 



(30) 



where k (k') is the wavevector of the incident (scattered) neutron, N is the total number 
of H 2 molecules in the target, x is the fraction of H 2 molecules which have odd J (i. e. 
are ortho molecules), and the subscripts indicate the partial cross sections due to ortho 
molecules, to ortho-para conversion, and to para-ortho conversion, respectively. When the 
sum over nuclear spin states is performed, these partial cross sections are given by 

4 J i= 0,Jf=l Z 

Si^, = \{b'f E P l KE-E % + E f )\{f\e^^{\n.p)\i)\ 2 
4 j t =i,j f =o A 

= i(60 a £ PAE - E t + E f )\(f\e iK ^ cos(i K ■ p)|.)| 2 , (31) 
z Ji=i,j f =x z 

where k, = k' — k, Rj is the center of mass of the jth H 2 molecule, and p is the vector 
displacement of one proton relative to the other proton in the H 2 molecule. To deal with 
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molecular orientations for molecules where J is at most unity, we write 

sir^K . p) = An 3l ^K P )Y.Y?(KYY?{p) , (32) 

where j n is a spherical Bessel function. We will assume that k is small enough that the term 
in j 2 can be neglected. Since it is not trivial to obtain a meaningful result which properly 
contains the Debye- Waller factor, we proceed simply, as follows. From the numerical solution 
on a mesh of points we obtain the family of wave functions (each one denoted \J,P;a)), 
for which J and P are good quantum numbers and a = 1, 2, 3, . . .. (In limiting cases one 
may replace the single index a by two indices N and 7, where N, the number of phonons, 
is nearly a good quantum number.) If we label the radial mesh points by k — 1, 2, 3..., then 
we have 



\J,P;a) 



1/2 

27rEi=_jE*|c^(fc^)l 2 r* 



(33) 



Here \r k ) (and later (r^|) is a wavefunction of unit amplitude at the position r k . Also the 
Cj*p(k,[/,ys are the set of coefficients (for fixed J, P, and a) which are obtained by the 
numerical solution of the (2 J + l)-component radial eigenvalue problem on a set of mesh 
points {rfc}. This discretized eigenvalue problem involves diagonalization of nonsymmetric 
matrix. (The radial equation gives rise to a Hermitian problem only if proper account is 
taken of the radial weight function.) The numerical program takes no account of any weight 
factor, but rather normalizes these wavefunctions by requiring that the sum of the squares 
of their coefficients be unity. Since we always wish to define inner products with a weight 
factor r-fc, we will explicitly include a factor r k when we take inner products. Then if X is a 
quantity which is local in r and r but may be off-diagonal in J and/or J z , we express its 
matrix element between such numerically obtained wavefunctions as 

(J',P';(3\X\J,P;a) 

jgjgr E W (A J z = fJ-'\X(r k , 4> r )\J, J z = /^>r fc cg;p,(fc > /x / )* C g,(fc,/x)e^- p+ ^-^ Q ^ 

= r TT72 ^ 

2tt Efc, M |cj Q p(M)| 2 niEfc>' \cj%'(k',n')\ 2 r k/ 

Here // assumes integer values between — J' and +J' and p integer values between — J and 
+J. 
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To evaluate the cross sections, the major problem is to evaluate the matrix element, 
which we may call < f\X\i >. We will not discuss all possible transitions (which are shown 
in Fig. 4 of I). Instead we will focus on the neutron energy loss spectrum due to a) para to 
ortho conversion and b) radial phonon creation on an ortho molecule. 



A. Para to Ortho Conversion 



For neutron energy loss due to para to ortho conversion, we need the matrix element 



X = 47ri 1 (^p)exp(^-R,)^^(/ t )^(p)^ 



(35) 



and for simplicity we consider the case when k is perpendicular to the cylindrical axis of 
the nanotube. Then we may as well place k along the local x axis. Also the center of mass 
of molecule j is at relative to the axis of the tube at position nf^ which contains the jth 



molecule. Then we may write 



(36) 



At low temperature the initial state (whose energy is denoted Ej) will be the ground state 
for J = and for a small value of P (which we denote Pi). Thus 



= \*[m\«P)?Z- X E e- E ^5[E -E i + E f ]M if , 
where Ef is the energy of the (J = 1) final state, and Z = Y^ i ex.p[—E i /(kT)] : and 



(37) 



M 



if 



It d( t>rEkJJ = l,Jz = »\YHP) ~ Yr\p)\J = 0,J Z = ^e^^^e^f-^+^ncjik^Yc^k) 



2vr 



Efc \ci(k)\ 2 r k Y, k4l \c f (k,fi)\ 2 r k 



1/2 



T l k[ J P f ~P t +i(. Kr k)cf(k, 1)* - J Pf -p.- 1 (Krk)cf(k,-l)*]r k Ci(k) 



^Hk\ci{k)\ 2 r k Y, k ^\cf{k^)\ 2 r k 

where J n {x) is a Bessel function and Ci(k) is the value of Cjp(fc) for the initial state and 
Cf(k) is the value of cffi P ,(k' , //) for the final state. 

We now discuss the qualitative meaning of this result. We calculated So^i for several 
temperatures and the results are plotted in Fig. 8a. At zero temperature the initial state 
has Pi — 0. So the energy loss is zero up to the cut-off energy which is the energy of para-to- 
ortho conversion. For small temperatures, the cross section does not appear discontinuously, 



(38) 
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but turns on rapidly over a range of energy of order kT. The first para-to-ortho (i.e. J = 
to J = 1) transition is observed at energies about 16.8 and 13.6 meV with approximately 
one-to-two intensity ratio, (corresponding to a splitting of 3.2 meV between J — 1,M = 
and J — 1, M — ±1 states, respectively). The center of gravity of the J = 1 levels gives the 
average para-ortho conversion energy, E c = 14.67 meV. This result represents only a small 
amount of downward shift of 0.03 meV from the free molecule value of 14.7 meV. 

We note that there are several neutron scattering experiments reporting the para-to-ortho 
transition P, 1(|. The observed splitting is about 1 meV, suggesting the idea that in those 
experiments hydrogen molecules were probably not inside the nanotubes. The calculated 
para-to-ortho splitting of 3.2 meV is slightly larger than the 2.6 meV splitting observed for 
H 2 on graphite 3] . 

In addition to the sharp para-to-ortho rotational transitions, Fig. 8a also indicates several 
broad radial phonon transitions at energies about 15 and 30 meV (similar values to those of 
H2 trapped in solid Ceo)- Finally we note that there are many lines in the spectrum due to 
transition between different tangential phonon states (i.e. P quantum number). However 
their observation could be problematic due to experimental energy resolution (which would 
be worse at high energies than FWHM of 0.5 meV used in Fig. 8). 



B. Ortho Cross Section 



Here we discuss the scattering from an ortho H 2 molecule. As before, the major problem 
is the calculation of the matrix element. In analogy with the previous results we write 

(39) 

where we neglect terms involving j2(|^p) and 

M if = |(/|e mrfcCOS ^|z>| 2 . 



Sx->! = -[b'jo^KptfZ-^e-^^SiE - E i + E f ]M if 

if 



Now for Mif we have in the notation of Eq. 



Mi 



if 



It d(j> r E fc ,„ e 1 ^ cos ^e^ p f- p ^r k c f (k, ^)*c t (k, /1) 



2tt 



Efc, M \ci(k, n) \ 2 r k J2 k ,a \cf(k, fi)\ 2 r k 



1/2 



E& )A , J Pf - Pl (^r k )r k c f (k, fi)*Ci(k, p) 



J2k,^ \ci(k, p)\ 2 r k Ek,n \cf(k,p)\ 2 r k 
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(40) 



(41) 



Figure 8b shows the calculated spectrum ortho cross section S\^i, indicating many transi- 
tions between a large number of states. It is possible to identify the radial phonon transitions 
for only one and two phonon states as indicated in the figure. On the other hand, the tran- 
sitions between tangential phonon states (i.e. states with different quantum number P) 
dominate the calculated spectrum, giving rise to many sharp peaks. Due to experimental 
energy resolution, it is probably not possible to observe the transition at high energies (say 
above 20 meV). However the resolution at energies below around 10 meV could be about 0.5 
meV (which is used in Fig. 8) and therefore it may be possible to observe these transitions. 

Finally we note that the tangential phonon transitions below 10 meV show a maximum 
near 3.6 to 4 meV. We can understand this by considering the condition that the neutron 
waveform become resonant with the wavefunction of H 2 molecule going around the circumfer- 
ence of the minimum of the Mexican hat. The phase change when the neutron passes through 
a diamater of the Maxican hat is 2kx , where r is the radius at which the Maxican hat po- 
tential is minimal. The phase change of the H 2 molecule going around half a circumference 
is 7rP. If we assume an initial state with Pi = 0, then the resonance condition is nPf = 2ktq. 
With k = 3A- 1 and r = 3.5A, we find P f « 6. Then E f - E { = h 2 P] ((2mr )~ 1 ) = 0.09P 2 
meV = 3.2 meV, in reasonable agreement with the numerical evaluation. 

VII. CONCLUSION 

We list the major conclusion from our study of H 2 molecules bound to nanotubes which 
we treat as smooth cylinders. 

• We have derived the analog of a radial equation for the Schrodinger equation for the 
translational and rotational motion of a molecule in cylindrical geometry. This formulation 
leads to classifying translation-rotation wavefunctions according to their properties under a 
global rotation of the molecule about the cylindrical axis. 

• Using this radial equation, the translation-rotation wavefunctions for a hydrogen 
molecule bound either inside or outside a nanotube can be obtained numerically. We also 
have developed simple toy models which quite accurately reproduce the numerical results, 
but have the advantage that they elucidate the nature of the translation-rotation dynamics. 

• Simple classical symmetry arguments fail to predict the correct degeneracies of 
translation-rotation wavefunctions. However, the quantum wavefunctions are easy to un- 
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derstand qualitatively. For instance, for a J = 1 molecule (such as ortho-H 2 ), one class of 
translation-rotation wavefunctions has the molecule in a J z = state (i. e. aligned along 
the axis of the nanotube) with no admixture from J z = ±1. This simplification is a result 
of the mirror plane perpendicular to axis of the cylinder. 

• We also suggest that neutron time-of-flight spectra could provide useful confirmation of 
our results. To that end we have calculated typical spectra that might be observed. These 
are shown in Fig. 8. 
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APPENDIX A: INCLUSION OF BOTH ANHARMONICITY AND 
TRANSLATION-ROTATION COUPLING 



Here we study the simultaneous effect of anharmonicity and translation rotation coupling 
for the two-phonon manifold within the toy model. From Table Ewe see that we may write 
the anharmonic Hamiltonian, Vah, which is independent of J z , as 

V AH = V 1+ ^7a 4 ^2|2,m = 0)(2,m = 0| - \2,m = 2)(2,m = 2\ - |2,m = -2)(2,m = -2|^A1) 

where Vo = 1527cr 4 /3 and |2,m) is a two-phonon wavefunction (with energy 3hu) as given 
in Table HI The translation rotation interaction may be written as 

V = a [3 J 2 -2] + V (A2) 

where 

V = ^(jl(x-iy) 2 + J 2 _(x + iy)^j . (A3) 
If we write the ground state wave function as 

|0) = N 0fl e-^ x2+y2)/a2 , (A4) 
then the two phonon eigenstates are 

|2, m = 0) = -(X 2 + Y 2 - 2)N m e-^ x2+Y2 "> (A5) 

and 

|2, m = ±2) = -^t=(X + iY) 2 N m e~^ x2+Y2) . (A6) 
2v 2 

where X = xjo and V = y/c. If ( ) indicates a ground state average, then we have the 
evaluation 

<0|(:r-^) 2 |2;m = 2> = a 2 [^-{X 2 + Y 2 - 2)(X - iY) 2 ^=(X + *F) 2 ^ 

= 4V2a 2 , (A7) 

so that 



V = 2v / 2>x 2 |j2;m = -2)(2;m = 0| + \2;m = 0)(2;m = 2|J J 2 + H. c. 



(A8) 



where H. c. indicates the Hermitian conjugate of the preceding term. We find the eigenvalues 
to be those of Table IIVI 
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APPENDIX B: CARTESIAN REPRESENTATION 

Here we rewrite the eigenf unctions in the Cartesian representation. Write 

1^ = 0) = \Z) , 

\Jz = l) = ~^\X + iY) , 

\J Z = -1) = ±=\X-iY) . 
Look at Table IIII1 Eigenfunctions with energy a ± 4<r 2 j3 are 

^ = ^{-^ x+iY ^ x - i y^^^ x ~ iY ^ x+i y^) 

So, apart from a phase factor, 

^+ = ^\xX + yY) 
i>- = ^=\xY-yX) . 



The eigenfunctions with energy a are 



1>± = -\\{x±iy){X±iY)) 



so that we may take the eigenfunctions to be 



= ±=\ x Y + yX) . 
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TABLE I: Two-Dimensional Harmonic Oscillator Wave functions 



E/{hu) 


dE / d~j) 1= o 






^(x,y) b 


1 


8a 4 





(0,0) 


1 


2 


24a 4 


1 




x + iy 


2 


24a 4 


-1 




x — iy 


3 


56a 4 





[ (2,0) + (0,2)] 


(r 2 /a 2 ) - 2 


3 


48a 4 


2 


[(2,0)-(0,2)+iV2(l,l)] 


(x + iy) 2 /a 2 


3 


48a 4 


-2 


[(2,0) - (0,2) -z^2(l,l)] 


(x - iyfja 2 


4 


96a 4 


1 


[V3(3,0) + (l,2)] 


[(x + iy)/a}[(r 2 a 2 )-A] 


4 


96a 4 


-1 


[V3(0,3) + (2,l)] 


[(x-iy)/a][(r 2 /a 2 )-4\ 


4 


80a 4 


3 


[(3, 0) - s/3(l, 2)] + iV3(2, 1) - i(0, 3) 


[(x + iy)/af 


4 


80a 4 


-3 


[(0, 3) - >/3(2, 1)] - z^3(2, 1) + i(0, 3) 


[{x - iy)/af 



(a) In cylindrical coordinates the ^-dependence is through the factor e tm ^ '. 

(b) The wave function contains, in addition to the factors listed, exp[— (x 2 + y 2 )/(4a 2 )] as well as 

a normalization factor. 
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TABLE II: Wavefunctions for Rotation- Vibration for J = 1, n x = n y = 0. 



Energy a 


pb 


Wavefunct 


ion 








= 0,n y 


= 






Jz = +l 


Jz=0 


Jz = -1 


-2a 








1 





a 


1 


1 








a 


-1 








1 



(a) We tabulate the energy relative to Two. 
(b) P defines the transformation of the wave function under a global rotation, as explained 

connection with Eq. 



TABLE III: Wavefunctions for Rotation- Vibration for J = 1, n x + n y = 1. 



Energy 3, 


pb 


Wavefunction 






n x 


= l,n y 


= 


n x 


= 0,n y 


= 1 






Jz = +l 


Jz =0 


Jz = "I 


J z = +l 


J 2 = 


J 2 = -1 


a + 4<7 2 





1 

2 





1 

2 


-¥ 





ii 


a - 4a 2 





1 
2 





1 

2 


-¥ 





-|i 


-2a 


1 





1 

V2 













-2a 


-1 





1 

V2 













a 


2 


1 

V2 
















a 


-2 








1 

V2 








"75* 



(a) We tabulate the energy relative to 2%lo. 
(b) P defines the transformation of the wave function under a global rotation, as explained 

connection with Eq. 
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TABLE IV: Wavefunctions for J = l,n x + n y = 2 with anharmonic (scaled with 7) and rotation- 
translation coupling (scaled with a). 



p 


Wavefunction a 


Energy b 







\ m = o ; M = 0) 




-2a + I67CJ 4 


1 


[|m = 


2;M = -1) + \m = 0; M = 


l>]A/2 


a + 127a 4 + AR 


-1 


[\m = 


-2;M = 1) + |m = 0; M = 


"1)]/V2 


a + 127a 4 + AR 


1 


[\m = 


2;M = -1) - \m = 0;M = 


1)]/V2 


a + 127a 4 - AR 


-1 


[\m = 


-2; M = 1) - \m = 0; M = 


-1)]/V2 


a + 12 7 cr 4 - AR 


2 




\ m = 2-M = 0) 




-2a + 870- 4 


-2 




|m = -2;M = 0) 




-2a + 87a 4 


3 




\ m = 2-M = 1) 




a + 87a 4 


-3 




\m = -2;M = -1) 




a + 87a 4 



a) Here M indicates a wavefunction for which J z = M. The states indicated by m are the phonon 

states in the cylindrical gauge and are listed in Tabled 

b) Here R = a/7 2 ct 8 + 2/3 2 er 4 . Also these are energies relative to 3Hlo + 4O70" , which is the energy 

for equally spaced levels extrapolating from the zero phonon and one phonon level. 
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TABLE V: Energy Levels in meV for the toy model for a ( J = 1) hydrogen molecule inside of a 
(9,0) tube, compared to numerical calculations based on Eq. (|lUj) : The parameters (in meV) are 
a = 2.82. at = 3.56, a 2 = 4.32, (3 1 cr 2 = 0.59, (3 2 a 2 = 0.55, 7a 4 = 0.29, and %uj = 27.36. 



p 


Toy Model Energy a ' b 


Energy a (Numeric) 











±1 


3a = 8.46 


8.46 





2a + a 1 + 4/3icr 2 + hu = 38.92 


39.01 





2a + ai- 4/?icj 2 + hto = 34.20 


34.29 


±1 


2a - 2ai + Huj = 25.88 


25.89 


±2 


2ao + «i + hu; = 36.56 


36.56 





2a - 2a 2 + 45 + 2Tiuj = 56.36 


55.67 


±1 


2a + a 2 + 3J + 4R + 2Tiuj = 71.68 


71.13 


±1 


2a + a 2 + 35 - AR + 2?iuj = 64.64 


64.12 


±2 


2a - 2a 2 + 87a 4 + 2huj = 54.04 


53.53 


±3 


2a + a 2 + 87a 4 + 2Tiu = 67.00 


66.35 



a) The zero of energy is taken to be the lowest P = level. 



b) Here R = ^ffa* + 2f3 2 a^. 
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FIG. 1: A schematic representation of a single wall carbon nanotube rope indicating three differ- 
ent absorption sites; namely (A) endohedral, (B) interstitial, and (C) external adsorption sites, 
respectively. 
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FIG. 2: Potential energy as a function of distance from tube center for a para hydrogen (i.e. J = 0) 
interacting with a (9,0) (top) and (10,10) (bottom) nanotube. The solid and dashed lines are for 
smooth and actual carbon nanotubes, respectively. For the actual carbon nanotubes, the value 
of z-axis is taken arbitrarily. The horizontal lines in the bottom panel indicate the radial phonon 
bound states in smooth nanotubes. The insets to the left panels are schematic plots of the parabolic 
(top) and Mexican hat potentials, respectively. 
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FIG. 3: (a) Potential energy as a hydrogen molecule is translated from the center of a (10,10) 
nanotube when H2 is oriented to be parallel to the tube axis (p), radially (r), and tangentially (t), 
respectively. Inset shows these configurations. The equilibrium distance and minimum potential 
depends strongly on the orientation of the H2 molecule, indicating strong rational-translational 
coupling, (b) The minimum potential energies, E r , E p , and Ef, for radial, parallel, and tangential 
orientation of H2 molecule for various nanotubes. 



31 




FIG. 4: Removal of the degeneracy in the energy level scheme of the one-phonon (J = 1) manifold 
according to the Hamiltonian of Eq. (|22|). The diagram labeled "SPHERE" is for a spherical 
molecule for which a = /3 = 0. That labeled "ROD" is for decoupled rotations and translations of 
a rod-like molecule for which That labeled "ROT-VIB" is for translation-rotation coupling 

with (3^0. 
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FIG. 5: Translation-rotation wavefunctions for a (J = 1) with one phonon when there is dynamical 
mixing of translations and rotations. Here the plane of the paper is the x — y plane and each 
figure eight represents an \X) or \Y) orientational wavefunction and the sign associated with each 
lobe of this p-like function is indicated. For the \Z) orientational function (which would have 
the figure eight coming out of the page) we indicate the sign of the lobe in front of the page. 
Each orientational wave function is multiplied by a translational wavefunction \x) or \y), where 
for instance \x) ~ x exp[— \{x/cj) 2 }. The presence of a phonon in the r a coordinate thus causes 
the wave function to be an odd function of v a , as one sees in the diagrams. One sees that the 
P = wavefunctions are invariant under rotation by tt/2. (In fact, they are angular invariants.) 
From the states labeled with nonzero values of P, one can form the complex linear combinations 
which transform like e lF ^ when the position and orientation of the molecule are simultaneously 
rotated about the symmetry axis. Although it is far from obvious, the two states labeled P = 2 
are degenerate in energy. The two P = states have different energy, in general. So quantum 
mechanics predicts the six state manifold to consist of two doublets (one for P = ±1 and one for 
P = ±2) and two singlets (for P = 0). 
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FIG. 6: Removal of the degeneracy in the energy level scheme of the two-phonon (J = 1) manifold 
according to the Hamiltonian of Eq. (|22|). The diagram labeled "SPHERE" is for a spherical 
molecule for which a = [3 = 0. That labeled "ROD" is for decoupled rotations and translations, 
but with a / 0. That labeled "ROT-VIB" is for anharmonic and translation-rotation coupling 
with (3^0. Here R = y/i*o* + 2/? 2 cr 4 . 
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FIG. 7: The energy levels of of an H2 molecule with (J = 0) (a) and (J = 1) (b) inside a 
(10, 10) nanotube versus quantum number P. The symbols and the dotted lines are obtained from 
numerics and the solid lines are fit based on the simple models as discussed in the text, indicating 
that a few of the lowest energy levels can be understood from these simple models. For the case 
of (J = 1) hydrogen (bottom), for each P we have now three energies, which are split by about 3 
meV (comparable to 2.6 meV observed for H2 on graphite). We note that for both (J = 0) and 
(J = 1) cases, only N = 0, 1, and 2 phonon levels can be safely identified. The inset to the lower 
panel gives the fitted values of the parameters of Eq. (|28[). In the lower panel, the energy is with 
respect to BJ(J + 1). 
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FIG. 8: The calculated neutron cross section (with arbitrary scaling) for para to ortho (a) and ortho 
to ortho (b) transitions at several temperatures. The neutron wavevector transfer k is taken to be 
perpendicular to the tube axis with magnitude of 3 A -1 . The peaks are broadened by Gaussians 
with FWHM of 0.5 meV. 
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